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ABSTRACT 

A  constitutive  two-phase  material  model  is  formulated,  which  satisfies  the  decoupling  of 
deviatoric  stress  and  pressure  response  required  by  many  hydrocodes.  Two-phase  materials 
are  assumed  to  be  non-dispersive.  Basis  of  the  present  model  is  a  decoupled  model  based  on  a 
general  Maxwell-type  viscoelastic  model.  The  proposed  formulation  provides  thermodynamic 
consistency  for  the  case  of  small  elastic  shear  strains  and  unrestricted  volumetric  response. 
The  model  is  verified  against  experimental  data  on  the  shock  wave  propagation  in  iron 
resulting  in  alpha-epsilon  phase  transition,  while  the  strength  effects  being  observed  at  the 
sameti  me. 
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Thermodynamically  Consistent  Decoupled  Strength 
Model  for  Two-Phase  Materials  with  Phase 

Transitions 


Executi  ve  Su  mmary 

In  order  to  increase  the  efficiency  of  modern  weapons  and  protection  systems,  novel 
technology  developments  deal  with  advanced  multi-phase  and  multi-component 
materials.  Among  them  are  filled  energetic  materials,  porous  and  multi-phase 
mitigants,  composite  materials,  including  nano-ceramics,  etc.  Behaviour  and  response 
of  a  variety  of  multi-phase  non-dispersive  materials,  including  porous  and  condensed 
phase-transforming  multi-phase  materials,  has  to  be  described  in  order  to  enhance  the 
technology.  This  poses  new  challenges  in  the  area  of  constitutive  multi-phase 
modelling  of  the  material  response  at  extreme  conditions.  Advanced  numerical  codes 
are  employed,  amongst  which  hydrocodes  are  particularly  suitable  because  they  may 
model  the  material  behaviour  in  a  wide  range  of  conditions  from  moderate  impact  and 
blast  loads  resulting  in  the  strength-sensitive  response  till  extreme  explosion  and 
hypervelocity  impact  loads  resulting  in  the  strength  negligible  response.  From  the 
physical  mechanism  of  deformation  viewpoint,  the  moderate  loads  are  associated  with 
the  el asto- plastic  deformation/  transition  and  the  extreme  loads  are  associated  with  the 
hydrodynamic  material  flows  at  extreme  pressures  and  temperatures.  The  elasto- 
plastic  response  in  the  defence  applications  is  needed  to  be  constitutive  (i.e.,  it  should 
take  into  account  the  rate  sensitivity  of  the  yield  stress)  and  requires  advanced  models. 
At  present,  the  multi-phase  material  models,  which  describe  both  the  strength  effects 
and  the  extreme  hydrodynamic  behaviour,  are  rare.  Therefore,  development  and 
implementation  of  advanced  models  in  commercially  available  hydrocodes  is  a  great 
challenge.  The  modelling  capability  in  DSTO  is  supported  by  two  major  hydrocodes, 
LSDYNA3D  and  CTH .  Hydrocode  LSDYNA3D  proved  its  efficiency  for  the  material 
response  to  moderate  deformations  when  a  good  resolution  of  contact  interfaces  is 
required.  However,  for  the  problems  such  as  counter-IED  (Improvised  Explosive 
Devices)  designs,  where  intensive  flows,  large  deformations,  and  multi-phase  mix-up 
are  common,  a  Eulerian  code  such  as  CTH  is  better  suited  for  this  model.  Both  codes 
require  material  models  decoupling  the  volumetric  and  shear  response,  which  are  not 
common  for  the  majority  of  thermodynamically  consistent  models  dealing  with  the 
wide  range  of  loads. 

Recently,  a  decoupled  model  that  describes  both  the  strength  effects  and 
hydrodynamic  behaviour  of  condensed  materials  has  been  developed  in  DSTO  [1]. 
However,  the  model  [1]  is  dealing  only  with  single  phase  materials.  The  present  work 
extends  the  model  to  the  case  of  two-phase  materials.  The  model  is  verified  against 
experimental  data  on  the  shock  response  of  iron  resulting  in  the  alpha-epsilon  phase 
transition.  At  the  same  time,  the  strength-level  response  is  observed  as  a  multi-step 
shock  wave  structure  with  the  steps  corresponding  to  both  the  alpha-epsilon  transition 


and  the  el  asto-plastic  transition  (elastic  precursor).  The  modelling  results  correlate  well 
with  the  experimental  results  available  in  literature. 
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1.  Introduction 


A  large  number  of  advanced  multi-phase  and  strength  material  models  have  been 
developed  in  the  last  decades.  Many  of  the  models  have  been  implemented  in  commercial 
hydrocodes  such  as  LS-DYNA  [1]  and  CTH  [2],  Many  of  the  strength  models  have  a 
decoupled  representation  separating  the  volumetric  (bulk)  response  managed  by 
hydrostatic  pressure  from  the  deviatoric  response  managed  by  the  shear  stresses.  A 
number  of  advanced  strength  models  that  are  attractive  from  their  thermodynamic 
consistency,  however,  are  not  decoupled  and  cannot  be  directly  used  for  the  hydrocode 
implementation.  In  order  to  keep  the  thermodynamic  correctness,  a  revision  of  one  of  the 
models  [3]  has  recently  been  conducted  and  a  decoupled  strength  model  has  been 
developed  [4]  and  implemented  in  the  CTH  hydrocode. 

The  present  work  generalises  the  model  [4]  to  the  case  of  two-phase  non-dispersive 
materials.  Such  materials  have  components/ phases  with  the  particle  velocities 
equilibrated  over  the  multi-phase  mixture.  Examples  of  these  materials  are  numerous  and 
include  porous  materials,  liquids  with  bubbles,  reaction  products  of  energetic  materials 
with  particles  small  enough  to  achieve  almost  instantaneous  equilibrium  by  velocity, 
phase-transforming  two-phase  condensed  materials,  etc.  The  present  extensions  of  the 
model  [4]  for  the  case  of  small  elastic  shear  deformations  and  unrestricted  bulk  response 
focuses  on  1)  formulation  of  the  model  into  the  form  that  is  in  agreement  with  the 
hydrocode  implementation  decoupling  convention;  and  2)  checking  out  the 
thermodynamic  correctness  of  the  model. 

Formulation  of  the  model  is  based  on  two  major  points:  i)  a  choice  of  thermodynamic 
variables  and  ii)  derivation  of  constitutive  equations  and  conservation  laws.  The  closing 
relations,  including  the  equation  of  state  (EOS),  are  obtained  naturally  from  the 
corresponding  relations  for  individual  phases.  While  implementing  the  model,  a  one¬ 
dimensional  in-house  wave  propagation  code  has  been  developed  employing  the 
Godunov  method  [5]  on  the  basis  of  the  Riemann  problem  procedure  [6]  for  the  bulk 
response  and  a  linearised  method  similar  to  the  method  of  characteristics  for  the  shear 
response. 

The  model  is  verified  against  data  available  in  literature  on  the  velocity  of  the  free  surface 
of  iron  target  in  the  plate  impact  experiments  at  using  the  shock  propagation  basic  set-up 
for  shock  compression  of  pure  iron.  The  calculations  demonstrate  both  alpha-epsilon  (bcc- 
to-hcp)  phase  transition  and  elastic  precursor  development  (the  elasto-plastic  transition), 
which  agrees  well  with  the  experiments. 


2.  Background  Constitutive  Model 

The  Maxwell-type  viscoelastic  model  [3]  reduced  to  the  case  of  small  shear  elastic 
deformations  is  used  in  the  following  form  [4] : 
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The  equations  are  written  down  in  the  Eulerian  coordinate  system  in  the  form  of  the 
balance  laws  with  external  exchange  terms.  The  source  terms  of  the  external  flows  mo,  n‘o, 
and  lo  are  zero  (turning  the  balance  laws  to  the  conservation  laws)  for  the  original  single 
phase  model  [4]  describing  a  closed  material  system.  We  will  denote  the  substantial 
derivative  d/dt  =  d/dt  +  uk  d/dxk.  The  prime  sign  will  refer  a  tensor  to  its  deviatoric  part,  so 
s ' j  =  Ey  -  Ett'Sy/3,  where  5,y  are  components  of  the  unit  tensor.  Here  e,y  =  s'y  / p  where  Ey  and 
Sij  are  components  of  the  tensor  of  small  elastic  deformations  and  deviatoric  stress  tensor, 
u,  are  velocity  components,  and  p  and  e  are  density  and  specific  internal  energy.  The 
velocity  module  u  is  determined  as  u2  =  w,:  •  Uj.  The  deviatoric  part  Sy  of  the  stress  tensor  a,y 
and  the  spherical  (negative)  part  p  of  the  stress  tensor  (pressure)  are  linked  as  follows 

Sij  =  Gy  +p  -8y  ,  p  =  -  akk/  3  .  (2) 

Pre-given  functions  x  (sy,  T)  and  e(p,  ejr  s),  where  T  and  s  are  temperature  and  specific 
entropy,  close  the  system  of  equations  of  the  model.  The  function  x  (s,y,  T)  is  obtained  from 
the  data  on  yield  limit  versus  strain  rate  (e.g.,  see  [7]  on  how  this  can  be  done)  and  the 
internal  energy  e(p,  ey,  s)  is  obtained  from  high-pressure  material  response  data  such  as 
Hugoniots  or  static  compression  data. 

The  thermodynamic  identity  (3)  for  the  present  model  (see  [4])  is: 

T  ds=  de  -  Sij  deji  +  p  dV,  (3) 

where  V  =  1/p  is  the  specific  volume.  Thus,  we  have  11  independent  thermodynamic 
variables,  namely,  e,y,  p,  and  s  and  11  dependent  variables:  s,y,  p,  and  T.  The  dependent 
variables  are  calculated  from  the  identity  (3)  as  follows: 

T=de/ds  ,  p  =  p2de/dp  ,  s,y  =de/deji  .  (4) 

From  the  symmetry  conditions  for  the  strain  and  stress  tensors  and  from  the  following 
deviatoric  conditions  for  tensor  deviators:  en  +  en  +  eyy  =  0  and  su  +  S22  +  S33  =  0,  the 
number  of  variables  to  be  calculated  is  reduced  to  10;  this  set  includes  the  thermodynamic 
variables  and  velocity  components. 
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3.  Thermodynamic  Variables 


For  a  two-phase  material  the  phases  are  characterised  by  density  p,  small  elastic  deviatoric 
deformations  ey,  specific  entropy  s,  plus  pressure  p,  deviatoric  stresses  sy,  and  temperature 
T;  the  second  half  of  the  variables  are  dependent  thermodynamic  variables.  Velocity  u  is 
the  same  for  both  phases  in  the  non-dispersive  materials.  Variables  for  the  corresponding 
phases  1  or  2  are  referred  by  superscripts  in  parentheses.  Both  phases  are  assumed  to  be 
compressible  and  to  have  individual  strength  and  other  material  properties.  The  volume 
and  mass  concentrations  of  the  first  phase  are  0  and  c;  they  are  defined  identically  to  the 
definitions  from  [6,  8-10]  with  similar  relations  for  averaged  density  p,  entropy  s,  pressure 
p,  and  internal  energy  e  : 


p  =  6  pF)  +  (1  -  0  )  pC 2) ,  p  =  p2ep  =  0  p( i)  +  (1  -  0  )  p( 2) , 

e  =  c  eF)  +  (1  -  c)  c(2)  ,  s  =  c  s(1)  +  (1  -  c)  s(2> . 


(5) 


For  deformations  we  assume  that  the  deformations  for  mixture  should  be  volume 
averaged  similarly  to  [11-13]: 


e'ij  =  0  eyi)  +  (1  -  0 )  e'y(2) 


Then,  using  the  relation  [4]  s'y  =  pey  and  the  definitions  [6,  8-10] 

p(1)  =  p  c/d  ,  p(2>  =  p(l  -  c)  / (1  -  0 ) ,  (6) 

we  have 


ey  =  c  etp  +  (1  -  c )  ef)  .  (7) 

Inter-phase  deformation  imbalance  is  introduced  as  follows:  p/\,;  =  c’,/1)  -  s',/2),  which 
results  in 


Xij  =  c  ey(  1)  /  0  -  (1  -  c)  ey(2)  /  (1  -  0  ) .  (8) 

The  entropy  disequilibrium  is  introduced  in  the  way  used  in  [6,  9]: 

X  =  c  sd)  -  (1  -  c)  s(2)  .  (9) 

The  definitions  (6)  and  (9)  enable  us  to  link  phase  entropies  and  densities  with 
corresponding  averaged  variables.  Following  these  definitions,  we  can  express  the 
independent  phase  state  variables  via  the  averaged  ones  as  follows: 

pF)  =  p  c/0 ,  p(2)  =  p(l  -  c)  / (1  -  0 ) ,  sF)  =  V2(s  +  x)  /c  ,  s(2)  =  Vi  (s  ~x)  / (1  -  c)  .  (10) 

Similarly,  the  link  between  e,/1),  e,/2),  and  ey  is 
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etP  =  6  (ej  +  Xy  (1-6))  /c,  efi)  =  (1-6)  (ey  -  Xy  6)  /  (1  -  c) .  (11) 

We  assume  that  the  internal  energy  dependencies,  EOSs,  for  each  phase  are  given  in  the 
form 

e(i)  =  e(i)(p(i),  eyO),  s(i))  ,  e(2)  =  e(2)(p(2),  e^),  s(2))  .  (12) 

Thus,  the  averaged  internal  energy  from  (5)  and  (12)  is  a  function  of  the  following 
parameters: 

e(p,  c,  6,  ej,  Xu,  s,x)  =  c  eO)  ( pc/ 6,  6(ey  +  A, ,(1-0  ))/c,  V2(s  +  j)/c)  + 

(13) 

+  (1  -  c)  eP)  (p(l-c)/(l-0 ),  (1  -  0)(c,-  XjjQ) / (1  -  c),  y2(s  -X)/(l  -  c))  , 

here  the  relations  (10)  and  (11)  have  been  used.  Summarising,  the  independent  variables 
that  fully  describe  the  strength  bearing  two-phase  material  are  p,  s,  u,  c,  0,  e(/,  Xy,  and  j.  For 
this  set  of  independent  variables  the  thermodynamic  identity  takes  the  following  form: 


T  ds=  de  +  p  dV  -  sy  de^  -  Qy  dXji  -  A  dc  - 1 1  d0  -  T1  rfy .  (14) 

Thermodynamic  identities  for  the  phases  can  be  used  for  calculation  of  averaged 
dependent  variables  from  (14)  as  follows: 

p  =  p~ep  =  p2(cem  +  (1  -  c)ei2)\  =0pm  +  (  1  -  0)p(2), 

=(«'"+(! -c)e'%  =es;;>+(i-e)p\  (is) 

T  =  e,=(cem  +  (\-c)em\  =cTm(sf,,\  +(l-e)r“(s121),  =  (r0)  +  Te,)/2 . 


Here,  the  chain  rule  and  definition  (13)  have  been  employed  when  differentiating  over 
averaged  independent  variables.  The  dependent  thermodynamic  variable  A  in  (14)  that  is 
responsible  for  the  change  of  phase  mass  is  traditionally  related  to  a  chemical  potential. 
Similarly  to  (15),  when  applying  derivatives  over  c,  X 0,  and  j  to  the  internal  energy  e  in 
(13),  we  have 


Cl)  (2) 

A  =  ec  -e("  +TiTT- r’V  -  4^  -  e*>  -  ?  +  T  “  >*<■»  +  e^sf  =  M,  -  fh , 


Q^e, 


.CD  _  „P) 


(16) 


n  =  =  -- 


p 


^-+k+0- 2fl)As](y»-sp) 


y  =  =(r(1)-r(2>)/2. 


Here  ///,  =  eW  +  /#)//#)  -  T(k')s<k'>  -  e^-Ws ,/h  is  the  Gibbs  energy  for  fc-th  phase  and,  thus,  the 
chemical  potential  A  associated  with  the  mass  concentration  change  is  proportional  to  the 
affinity  of  the  Gibbs  energy.  The  potentials  Qy  and  VE  associated  with  the  variations  of  the 
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strain  imbalance  A ,y  and  the  entropy  disequilibrium  j  are  proportional  to  the  affinities  of 
the  deviatoric  stresses  s;y  and  temperature  T.  From  (16),  the  potential  n  associated  with  the 
variation  of  the  volume  concentration  can  be  rewritten  as  follows 


,0) 


,(2) 


n  =  ee  =-- 


,(!) 


y 


(. 


,(2)  )_42*( P  e' 


(2)„(2) 
y 


-;-rsV  -n"-;-rsy 


This  means  that  the  corresponding  affinity  is  linked  with  the  pressure  affinity  and  with  the 
affinities  for  the  work  of  elastic  shear  strains  and  that  for  the  strain  (the  last  affinity  is 
simply  the  strain  imbalance).  Using  the  definition  of  the  strain  imbalance,  the  last 
expression  can  be  rewritten  as  follows 

^(1)  -  (  p(i)41}  )  (4I}  -  ^ (2) + ( >^(2)42) )  (42)  -  ^ ) 

n  =  Cg  -  • 

p 


Thus,  the  volume  change  associated  potential  IT  is  proportional  to  the  pressure  affinity 
and  the  affinity  of  (pMe^-^M  -  s,y)  that  is  the  shear  strain  work  of  the  fluctuations  of  the 
deviatoric  stresses  over  the  averaged  stresses. 


4.  Two-Phase  Constitutive  Model 


In  order  to  derive  differential  equations  of  the  model,  we  will  employ  the  equations  (1)  of 
the  background  model  to  the  individual  phases.  In  doing  so,  the  source  (exchange)  terms 
that  are  denoted  for  the  first  phase  by  mo,  n\  and  lo  will  be  replaced  by  -mo,  -Wo,  and  -Iq  for 
the  second  phase.  This  provides  the  balance  of  mass,  momentum,  and  energy  in  the  closed 
two-phase  material  system.  For  a  representative  volume  containing  several  phases,  the 
density  and  pressure  are  characterised  by  partial  densities  pF) p  =  p-c  =  0-pF)  and  p(2) p  =  p-(l- 
c)  =  (l-0)-p(2),  and  by  partial  pressures  pF)p  =  0-pF)  and  pF) p  =  (l-0)-p(2)  (see  [9]).  The 
pressures  are  obtained  from  (4)  for  each  phase,  where  partial  densities  are  replacing  the 
actual  densities  pF)  and  p(2)  in  EOS.  Similarly,  partial  characteristics  for  deformations  and 
stresses  can  be  introduced,  resulting  in  e,yF)p  =  e,f)/  0,  e,y(2)p  =  e,y(2)/ (1  -  0  )  ,  and  s,yF)p  =  0-SyF), 

Sy(2)p  =  (1  -  0  )  -Sif). 

Summing  up  the  mass,  momentum,  and  energy  balance  laws  (1)  for  the  first  and  second 
phase,  we  have 
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dp  |  dpuk  _  Q  dpui  |  ojp11,",  -sy)  |  op  _Q 
dt  dxk  ’  dt  dxj  dxt 

dp[e  +  u2/2  )  d  puj  [e  +  )- 5 .. m,  +  pu  .  _  _ 

Sxy. 


(17) 


For  the  volume  concentration  6  the  constitutive  equation  is  taken  in  the  following 
traditional  form 


dp  6  dpOiij 

dt  dxj 


~PW  ■ 


(18) 


Equation  (18)  can  also  be  rewritten  as  Ad /dt  =  -  ip  with  the  use  of  the  mass  conservation 
law  from  (17). 


Using  mass  balance  equation  for  the  partial  density  p(9 p  =  p-c  of  the  first  phase,  we  can 
associate  the  rate  of  change  of  the  mass  concentration  parameter  c  with  a  rate  function  cp. 
Thus,  the  continuity  equation  for  p(1)p  gives: 


dpc  dpcuj 
dt  dxj 


pep. 


(19) 


From  the  constitutive  equation  for  c  of  (1),  it  is  obvious  that  mo  =  -  pep.  Using  equations  (19) 
and  (18),  we  can  derive  auxiliary  equations  for  pO)  and  p(2): 


dp (1)  dp^Uj  _  m0  +  p(l V  dp{2)  dp(  \ij  _  m0  +  p(2)y/ 

dt  dxj  9  ’  dt  dxj  1  -  6 


(20) 


The  equations  (20)  are  needed  for  derivation  of  constitutive  equations  for  strains  and  strain 
imbalances.  In  particular,  the  strain  constitutive  equations  from  (1)  for  the  phases  are 


dp(x)pefp  dp[i)pukefp  dU/2  dl 


dxk 

dxj 

dxt 

dui/ 

duj/ 
°  /2 

dUk/ 

f  /3 

dxk 

dxj 

<9x; 

Sx, 

p^pefp 


Keeping  in  mind  that  p(*)pei/0)p  =  p(%y(  1)  and  p(2)PCy(2)p  =  p(2)ey(2),  and  using  (20),  from  (21)  we 
have: 
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m0  +  /7('V^ 


e 


+p 


fi  U 


5  */ 


y 


dt 


fir, 


fir, 


-  +  ■ 


dx. 


-A=- 


UiUO 


-  e 


(2) 


m0  +  p(2)Y^ 


i-e 


+p 


/ 


d  J/ 


(22) 


fiM/ 


y 


dt 


dx , 


fir. 


■  +  ■ 


fir 


7  ^(2)^(2) 

3_.<y..=_^_SL  =  _^(2) 

11  T  11 


Now,  we  can  reduce  (22)  if  multiplying  the  equations  from  (22)  by  pc/ p(9  and  p(l-c)/p(2), 
respectively,  and  summing  them  up  separately  with  equation  (19)  for  c  (with  mo  replacing 
the  right  hand  side)  multiplied  by  c,/1)  and  with  a  similar  equation  for  1  -  c  multiplied  by 
e,/2),  respectively.  Using  (6),  the  resulting  equations  are 


dpcef)  dpce^uk 


dt 


dx,. 


-e 


r  dU'L  fiV  fiM*/ 

72  72  + — 73  •  £ 


fir. 


fir,  fir. 


/>c4}  a)  pi 

- py)  4V 


fip(l-c)ej2)  fip(l  -c)42,aU 


dt 

-M) 


dxk 

r  du‘/,  dUj/i  dUk/ 

— — — +^1.^ 


(23) 


fir, 


fir,  fir. 


p(l-c)« 


j(2) 

—  +  p(2)efy/  . 


Summing  up  equations  of  (23)  and  using  (7),  we  have  a  constitutive  equation  for  ey : 


dpej  |  dpe,jUk 

dt  dxk 


•  Sy  =  " +  0-  0 W ’)  -  P7  > 


(24) 


here  the  notations  from  (22):  cpyiX)  =  p  [e,y  +  (1  -  0)  A//]/ n  and  cpy(2)  =  p  [ey  -  6  Xy]/ T2  are  used. 
Direct  subtraction  of  the  second  equation  from  the  first  one  in  (21)  gives  the  following 
constitutive  equation  for  the  parameter  Xy  defined  in  (8) 


dpXy  |  dpXpp 
dt  dxk 


(25) 


Derivation  of  a  constitutive  equation  for  j  is  conducted  similarly  to  [9].  Specifically,  we  use 
the  energy  equations  for  phases  that  take  the  following  form  after  reductions  with  mass 
and  momentum  balance  laws 
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P 


(' )p  —  +  P{l)p  ^T~ ~ 4)p  ~r±  =  Io- Uk”o  - m0 (g(1)  - u/2 


p 


(2  )p 


dt 

de(2) 

dt 


du , 

l  P 

'J  dx, 

du  . 

ip 

',J  dx.. 


dxk 

+  pd-)p^_  s(2 )P  Ad  =  -/0  +  M  +  m  ( e 
dx. 


(2) 


2/ 


(26) 


Using  the  thermodynamic  identities  for  each  phase  written  in  the  actual  (non-partial) 
variables,  the  mass  balance  laws  (20),  the  momentum  balance  laws  for  the  phases  and 
constitutive  equations  (22),  from  (26)  we  can  find  rates  for  phase  entropies: 


p(  1)7(1)  aS  _  _ 


(  (1)  ,  /  (1)  (1)^1 
e  +  Zn®-S»e« 

V  /  P 


dt 

P[\ 


_ />); 

e 


(2)71(2) 


p^T 


ds(2)  l0+m0U/-m0 


6 


'e^  +  py, 

(  +  A  (2)  Sij  eij 

v  /  P  J 


■  + 


-«*«  0 


dt 

P[1)Y 


r 


y  y  / 

/  P  / 


1-0 


p(2)42)42) 


1-0 


(27) 


Multiplying  the  rate  equations  (27)  by  0/T  d)  and  (1-0)/ T  (2)  and  summing  them  up  with 
the  mass  balance  equations  for  phases  multiplied  by  sd)  and  s(2),  we  have 


g/AA  a(1A(A, 

df  fit* 


(28) 


Here  /r  V  =  /(/.  -  m2/  2  where  /4  is  the  Gibbs  energy  defined  in  the  previous  section. 

Summing  up  and  subtracting  the  equations  in  (28),  we  can  obtain  the  entropy  balance  law 
for  the  averaged  entropy  defined  in  (5) 
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dps  dpsuj 
dt  dxj 


-  R{+  R2 , 


(29) 


and  the  following  constitutive  equation  for  the  entropy  disequilibrium  parameter  j 
defined  in  (9) 


dpx  |  dPXUj 

dt  dxj 


=  Rx 


-r2 


—pco  . 


(30) 


Thus,  the  system  (17)  and  the  constitutive  equations  (18),  (19),  (24),  (25),  and  (30)  generate 
a  system  describing  behaviour  of  a  two-phase  non-dispersive  strength  bearing  material. 
The  conservation  laws  are  closed  with  the  equations  of  state  (12)  for  each  phase,  resulting 
in  EOS  (13)  for  the  mixture,  where  parameters  of  the  phases  are  connected  with  the 
averaged  parameters  by  relations  (10)  and  (11).  The  constitutive  equations  are  closed  with 
selection  of  the  constitutive  rates  (the  functions  cp,  ip,  and  o),  which  will  be  discussed  and 
specified  in  subsequent  sections. 


5.  Constitutive  Relations  and  EOS 


Similarly  to  [4],  we  choose  in  the  present  work  a  brief  form  of  a  Mie-Gruneisen-type  EOS 
suggested  in  [14]  and  modified  to  the  set  of  variables  in  the  decoupled  representation: 


,ey,s 


) 


2  a~n 


■(<?“•  -l)2  +  2plb;Sf’*1d  +c,T„S” 


1 


(31) 


This  EOS  is  used  for  each  individual  phase  of  the  two-phase  material.  Here  ao,  (3o,  and  yo 
(Gruneisen  coefficient)  are  material  constants,  cv  is  the  thermal  capacity,  po  is  initial 
density,  and  5  =  p/po.  It  should  be  noted  that  the  EOS  is  specific  to  individual  phase  with 
the  phase-specific  constants.  The  independent  variables  in  (31)  refer  to  /.)<'),  e,/1 ),  and  si1)  for 
the  first  phase  and  to  p(2),  e,/2),  and  s(2)  for  the  second  phase.  The  constant  ao  is  the  bulk 
sound  velocity  that  is  linked  with  the  longitudinal  and  shear  sound  velocities  Co  and  bo  as 
follows 


flo2  =  Co2 -Abo2/ 3  , 


and  d  is  the  second  invariant  of  the  strain  deviator: 

d  —  6ij  Zijl'2.  . 


Shear  stresses  are  calculated  from  (31)  by  the  rule  (15)  as 
sy  =  eey  =eddeij  =edeij  =  2pX^+2ev  =2pbi8Popej 


follows 
=  2ppeij 


(32) 


(33) 
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where  p  =  p-b2  =  p  ■ba25\'(>  is  the  shear  modulus  associated  with  the  shear  sound  velocity  b. 

Similarly  to  [4],  the  relaxation  time  functions  n  and  T2  in  the  constitutive  equations  for  the 
strains  e,/1)  and  e,/2>  are  chosen  in  the  unified  form  [15]: 


t(p,s,T)=  t0 


exp| 


(D0+Hel 


N0  +  Ms 


(34) 


Here  To,  Do,  H,  No,  M  are  material  constants  that  are  specific  to  individual  phases  as  well. 
The  stress  invariant  (s')2  =  s;y  -Sij  and  the  variable  e  is  associated  with  the  elastic  portion  of 
the  deformation  e  =  s'/ (2p).  The  material  constants  are  found  with  the  routine  procedure 
[7]  employing  experimental  data  on  yield  stress  versus  strain  rate. 


6.  Thermodynamic  Correctness  of  The  Model 

Verification  of  thermodynamic  correctness  we  start  with  expansion  of  the  energy 
conservation  law  in  (17)  reduced  with  the  mass  conservation  law  to  the  following 


de  du;  du dsH  du,  dp 
p —  +  pu .  — L  -  s,,  — L  -  u  — L  +  p — -  +  u,  -2—  =  0  . 
dt  dt  J  dxj  dxj  dxk  dxk 


(35) 


Using  the  mass  and  momentum  conservation  laws  from  (17)  and  the  chain  rule  applied  to 
(13),  we  transform  the  energy  equation  (35)  as  follows 


dp 

dt 


dc 

dt 


Pep  —  +  Pec  —  +  Peo  —  +  Pee, 


cW_ 

dt 


de 

ij 

dt 


+  Pex„ 


d/ 1 , 


ds 


dx 


— -  +  pe s  —  +  pev  — - 
dt  dt  x  dt 


du,  duk  A 

-s„ — p — -  =  0  . 
J  dxj  dxk 


(36) 


Using  again  the  mass  conservation  law  (17)  and  the  constitutive  equations  (19),  (18),  (24), 
(25),  and  (30)  for  c,  0,  etj,  Xij,  and  y,  the  equation  (36)  is  reduced  to  the  following 


Am0  -  pTlip  -  (dcp/1)  +  (1  -  d)cpiP)sij  -  pXijSijxp  -  ( cpi -  cpi/2))Qij  + 
+  pT-ds/dt  +  (R\  -  R2)(T( i)  -  T(2))/2  =  0  . 


(37) 


In  order  to  check  out  consistency  of  the  derivation  we  replace  the  particle  derivative  of 
entropy  by  equation  from  the  balance  law  (29)  and  the  potentials  A,  n,  and  Q,,  by  their 
expressions  from  (16): 


mo(/D)  -  /d2))  +  (pCO  -  p( 2)  -  p(1)e!;(1)Sy(1)  +  p(2)ey(2)Sy(2)  +  pXijSij)xp 
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-  ( dcpp  +  (1  -  9)cpij(2))sij  -  pXijSijxp  -  ( (pif !)  -  cpp)6(  1  -  0){siP  -  Sif2))  +  (38) 

+  (Ri  -  Ri)  (TO)  -  T( 2))/2  +  (Ri  +  R2)  (TO)  +  TO))/ 2  =  0  , 

here  the  corollary  of  the  definition  for  n  from  subsection  3  has  been  used.  Expansion  of  the 
last  two  terms  in  (38)  gives 

RiTO)  +  R2T(  2)  =  -  mo(jiO)  -  p))  +  (pO)  -  p0)  +  pO)ef)sf)  -  p(2kpsp)ip  + 


+  ecpijO)SijO)  +  (1  -  %y(2)Si/2)  . 


Replacing  the  last  two  terms  in  (38)  by  this  expansion,  we  obtain  an  identity. 

In  order  to  obtain  the  entropy  dissipation  term,  we  calculate  the  particle  derivative  of 
entropy  from  (37): 


pT ^  =  P<PA  +  +  (6(pf  + 11  “  #W2)K  +  PX  ‘js uV  + 

+ [vf  -  <pf ]  M1  -  °  M ]  -  sf )  ■ +  p°^  • 


Expanding  Sy  from  (15)  and  n  from  the  corollary  used  above,  we  have 


ds 

pT—  =  pcpA  +  py/U0  - pAijsijy/  +  pco^1  +pAijsijy/  + 


[o<pf  +  ('-  eWte4]  +  •'”  °)sf)+ ~  °tsf  ~ sf)  • 


(39) 


Transformation  of  the  last  two  terms  gives  dcppsp  +  (1  -  9)cppsp.  Thus,  the  equation 
(39)  takes  the  following  form 


PT~p  ~  PVK  +  P ^no  +  P^  +  e<pfsf  +  &  “  e)<pt\ 


(2) 

ij  ~ij 


(40) 


It  is  easy  see  that  the  products  cppsp  and  cppsp  are  proportional  to  epsp  and  epsp 
with  positive  coefficients.  These  terms  are  normally  positive  because  the  stress  deviators 
are  proportional  to  the  strain  deviators  with  a  coefficient  called  the  shear  modulus  as  in 
(33).  Thus,  in  order  to  achieve  non-negative  dissipation,  the  right-hand  side  of  (40)  can  be 
made  non-negative  with  the  following  choice  of  the  constitutive  functions  cp,  ip,  and  ox 

cp  =  A-cpo  ,  tp  =  Uo-ipo  ,  co  =  T'-fflo  ,  (41) 

where  cpo,  ipo,  and  coo  are  arbitrary  nonnegative  functions. 

One  case  of  practical  significance  is  the  one-dimensional  plane  set-up.  For  this  case,  we 
check  out  correctness  separately.  In  this  case,  the  system  (17),  (18),  (19),  (24),  (25),  (30)  takes 
the  following  form: 
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S_P  +  Spu_  =  ^ 
dt  dx 


dpu  d( 


pu  -  sr 


dt 


dx 


dx 


dp\e  + 


pu\e  +  u  /) )-  s1  u  +  pu 


dt 

dpO  dpOu 

+ - =  -PY  , 


-  =  0 


dt 


dx 


dx 

dpc  dpcu 

■  + - =  ~pq>  , 


dt 


dx 


(42) 


dt  dx  3  dx 
dpi  j  dpi  jW 


dt 


■+"aT=_Wll-"w)' 


Spy  Spy 
f  X  -  +  —  =  -pa>  . 


dt 


dx 


Here,  the  tensor  component  associated  with  the  x-direction  is  indexed  with  'T  and  index 
for  vector  variables  associated  with  the  same  direction  is  omitted.  Similarly  to  (36),  using 
(42)  for  expansion  of  the  energy  equation  we  have 


dc  dG  de,  dl ,  ds  dy  du 

per  —  +  peft  —  +  pep  — L  +  pe,  — L  +  pe  —  +  pe  —  -  s,  —  =  0  . 
dt  dt  1  dt  1  dt  dt  dt  dx 


This  equation  is  reduced  with  the  constitutive  equations  from  (42)  as  follows 


-ec<p-egY  +  - 


2  du 


-eza>- 


P 
sx  du 
p  dx 


,  .  (0(Px] +  ^-0)(p\2))~  pl,Y  [<Pi]  ~(p12))+T 

3  dx  J  p  dt 

=  0  . 


ds 


(43) 


The  equation  (43)  gives 


s,  =  —  e„ 


eccp  +egy/  [(%(1 1  +  (l  -  0)(p\1] )  +  pi  xy/\ +  ^ —  (pp  -  ^(2) ) 

cis  p  p 

dt  T 


.(44) 


It  is  easy  to  see  that  due  to  the  tensor  symmetry  restrictions  resulting  in  d  =  3ei2/ 4,  the 
shear  stress  defined  in  (44)  is  reduced  to  the  same  equation  (33)  in  the  form  si  =  2ppei.  The 
entropy  balance  equation  (44)  is  similar  to  the  entropy  balance  equation  (40)  in  the  general 
case,  when  taking  into  account  the  symmetry  induced  relations  s,,e,;  =  3si  t’i  /  2  and  Su(p,,  = 
3si(pi/ 2  that  are  valid  for  both  the  averaged  and  phase  variables. 

For  checking  out  hyperbolicity  of  the  system,  it  must  be  written  down  in  the  characteristic 
form.  The  present  system  in  the  corresponding  form  is  obtained  below  from  the  system 
(17),  (18),  (19),  (24),  (25),  and  (30)  as  follows 
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dp  dll, 

—  +  p — - 
dt  dxk 

deg  1 
dt  2  p 

dc 


=  0, 


dut  1  dsy  1  dp 


dt 


dt 


=  ~<P, 


8ui 

dxj 

dd 

dt 


du ,  2  5m 


- L  +  —^L  =  o, 

p  dxj  p  dxt 

\  t ^  (i) 


fix,  3  fix. 


■8„ 


K0+ (i-»K(2>) 


P 


(45) 


dAj  _ 

dt 


(<Py]  -  <pf] 


) 


P 


dX_ 

dt 


=  . 


+  i?2 


dt 


P 


It  is  obvious  that  structure  of  the  characteristic  system  (45)  is  coincident  with  the  structure 
for  the  single-phase  background  model  (1)  with  several  additional  contact  characteristics 
for  eigenvectors  corresponding  to  the  trajectory  evolution  of  variables  c,  9,  Xg,  and  j. 

Thus,  the  hyperbolicity  conditions  for  the  present  model  are  identical  to  those  from  [4] : 


4  ——z~ (s j )  =  a2  +  —  b2  =  c  >  0  . 
3 p1  V  3 


(46) 


However,  pressure  and  shear  stress  in  (46)  are  calculated  for  the  two-phase  mixture  using 
(15)  and  (13).  Nevertheless,  because  of  the  thermodynamical  correctness  of  the  model  the 
sound  velocities  are  calculated  with  the  mixture  rules  as  described,  e.g.,  in  [6,  10], 
Therefore,  hyperbolicity  of  the  system  is  achieved  with  the  traditional  requirement  of 
positiveness  for  the  expression,  representing  square  of  the  longitudinal  sound  velocity  c2. 


7.  Specification  of  Constitutive  Relations 

This  report  is  dealing  with  the  phase  transformation  study  for  condensed  two-phase 
materials.  An  essential  feature  of  this  case  is  that  both  phases  are  condensed  if  comparing 
with  the  porous  materials  analysed  earlier  when  strength  is  neglected  [6,  9],  Another 
feature,  which  differs  the  present  case  from  the  porous  material  analysis,  is  the  combined 
evolution  of  the  mass  concentration  due  to  the  mass  exchange  driven  by  the  phase 
transition  and  of  the  volume  concentration  due  to  different  compressibilities  of  the  phases. 
For  the  porous  material,  the  mass  exchange  did  not  occur,  which  was  tabulated  by  the 
choice  (p  =  0,  therefore,  the  mass  concentration  did  not  change  in  a  material  volume  [6,  9] . 

Because  of  this  two-fold  volume  variation,  we  need  to  take  into  account  the  volume 
exchange  due  to  possible  phase  transition.  Let  us  consider,  at  a  moment  of  time,  a 
representative  volume  V  of  mass  m  for  the  two-phase  mixture,  which  contains  the  both 
phases  with  a  volume  VW  of  the  first  phase  material  and  with  a  volume  V<2>  of  the  second 
phase  material.  Corresponding  masses  are  md)  and  m12).  In  order  to  calculate  the  change  of 
volume  and  mass  during  the  phase  transition  process,  we  need  to  rewrite  the  constitutive 
equations  for  mass  and  volume  concentrations  (18)  and  (19)  in  terms  of  m  and  V  as  follows 
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dyM 

dt 


-Vy/ 


dm ^ 

- =  —m  (p 

dt 


(47) 


Then,  the  volume  change  AW)  of  the  first  phase  within  a  small  time  increment  At  is 
composed  of  the  following  two  contributions: 

AW)  =  A  +  AVd1) , 

where  the  mass  change  contribution  AW(1)  =  Ami1) /pi1)  is  calculated  from  the  second 
equation  of  (47): 


a,  AT/W  -  _ -  m(PAt  -  y(l)(PAt 


Am  =  m  At  =  -mcp  At  ,  A  Vr/l 


P 


«  -  B%( 0 


and  the  volume  change  contribution  AW1)  is  calculated  from  the  first  equation  of  (47): 
AW1)  =  -  VipAt  . 

Thus,  the  total  volume  change  is 


AF{,)  =-Vy/ At - 


V(l)(pAt  1  AV 


(i) 


V  At 


cpd 

-V - =  ~¥, 

c 


Therefore,  for  a  prescribed  volume-change  kinetic  ip  driven  by  the  description  of  (18),  the 
constitutive  function  ipt  associated  with  the  total  volume  change  is  determined  as  follows 

ipt=  ip  +  dcp  /  c  , 

here  ip  =  TTo  ipo,  according  to  (41). 


Further  on,  we  specify  the  non-negative  function  cpo,  ipo,  and  oo  from  (41).  The  function  Mo 
can  be  selected  similarly  to  [6,  9]: 


_0  h  T(l)+T(2) 

a>0  p 


h  —  As 


k20  +  k]  (l  -  o) 


(48) 


Here  the  inter-phase  heat  transfer  coefficient  h  is  determined  according  to  [6,  9],  d  is  a 
characteristic  dimension  (typical  grain  size  for  a  phase).  As  is  a  dimensionless  factor 
(selected  as  a  value  of  the  same  order  as  previously.  As  =  300),  and  k\,  fe  are  the  thermal 
conductivities  of  the  phases.  The  function  ipo  is  chosen  in  the  same  form  as  [6,  9]: 

xp0  =  e(i-e)e0(p)  .  (49) 

The  multiplier  0(1  -  0)  keeps  the  parameter  0  within  the  [0,  1]  range.  A  function  9o(p ) 
responsible  for  the  compression  response  is  taken  as 
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0O(p)=Ao 


f  pT 


VPoJ 


For  the  mass  transformation  rate  cpo  we  take  the  Arrhenius  kinetic: 


(po  =B0 


VCmax  Ce  J 


-exp 
P  0 


Ugct-YsPoP 

RT 


(50) 


The  multiplier  with  parameter  c  keeps  this  parameter  positive  with  a  small  gap  c£:  c  >  ce. 
The  constants  R  is  the  gas  constant  and  po  is  the  averaged  reference  density. 

These  functions  cpo,  ipo,  and  Oo  from  (48-50)  will  be  used  in  the  subsequent  subsection  for 
the  shock  wave  propagation  calculations. 


8.  Phase  Transition  Calculations 


In  the  present  subsection  we  consider  an  example  of  phase  transition  in  pure  Iron  (Armco 
Iron)  loaded  by  shock  wave.  This  transition,  which  is  usually  observed  statically  at 
pressure  above  11-13  GPa,  is  associated  crystallographically  with  the  bcc-to-hcp  transition 
of  the  metal. 


Table  1:  Constants  for  the  constitutive  relations  and  EOS 


EOS  constants 

p0/  g/cm3 

ao, 

km/  s 

bo,  km/ s 

ao 

fo 

r° 

°Vr  J/g/K 

a-iron 

7.84 

4.637 

2.866 

1.73 

2.7 

1.735 

0.4415 

e-iron 

8.3 

4.67 

3. 

1.7 

2.7 

2.0 

0.4459 

Constitutive  Equation  constants 

h, 

W/m/K 

hi, 

W/m/K 

Ao, 

g/J/s 

no 

Bo, 

g/J/s 

mo 

d, 

pm 

Ce 

Cmax 

Uact/ 

KJ/mol 

r» 

g/mol 

80 

50 

1 

1.05 

3 

1.5 

10 

0.02 

0.9 

26.77 

0.127 

Elasto-plastic  transition  constants 

dei/dt,  s-1 

de2/dt,  s-1 

Yi,  GPa 

Y2,  GPa 

a-iron 

lO-3 

105 

0.2 

0.8 

e-iron 

lO-3 

105 

0.8 

3.2 

7o,  S 

Do/  GPa. 

H, 

GPa 

No 

M 

a-iron 

0.14 

4.174 

0 

106 

1011 

e-iron 

1.993 

16.7 

0 

106 

ion 
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To  specify  EOSs  for  (12)  in  the  form  of  (31)  for  a-  and  e-phases,  we  fit  constants  in  the 
EOSs  to  the  data  on  a-  and  e-iron  published  in  [16,  17].  The  results  are  summarised  in 
Table  1.  The  constants  chosen  for  the  constitutive  functions  cpo,  ipo,  and  coo  of  (48-50)  are 
stated  in  the  same  Table. 

The  constants  fitted  for  the  elastic-plastic  transition  kinetic  (34)  are  calculated  from  data 
taken  for  a-iron  from  [18]  and  stated  as  two  yield  limits  Yi  and  Y2  against  two  strain  rates. 
The  strength  data  for  e-iron  are  not  common;  indication  of  the  strength  increase  till  3  GPa 
reported  in  [19]  have  been  used  in  the  Table.  However,  for  the  majority  of  the  calculations, 
the  same  strength  values  have  been  used  for  the  both  phases  except  for  specially  stipulated 
cases.  Initial  density  of  the  iron  samples  in  the  calculations  is  chosen  to  be  po  =  7.86  g/ cm3. 


Figure  1:  Evolution  of  shock  profiles  at  lZo=0.1  km/s 

First,  the  test  calculations  have  been  conducted  for  shock  wave  propagation  through  a 
2  cm-sample  for  a  given  velocity  Uo  behind  the  shock  front.  This  set-up  may  be  seen  as  an 
analysis  of  the  target's  behaviour  in  the  plate  collision  one-dimensional  set-up  between 
two  iron  plates  with  the  flyer  plate  impact  velocity  of  2lfo.  The  variable  profile  in  these  test 
calculations  within  the  sample  thickness  are  drawn  at  equal  intervals  up  to  a  pre-selected 
time  Tf.  The  profiles  are  shown  for  velocity  u,  the  longitudinal  stress  Ox,  the  mass 
concentration  of  the  a-phase,  and  temperature  T. 

Results  of  the  first  calculation  at  Uo  =  0.1  km/ s  are  shown  in  Fig.  1.  The  final  drawn  profile 
6  in  Fig.  1  corresponds  to  Tf  =  3.1  psec.  The  calculation  demonstrates  that  the  elasto-plastic 
transition  occurs  in  this  case,  whereas  temperature  varies  insignificantly  from  its  initial 
value  To  =  293°K  and  the  mass  concentration  does  not  vary  at  all,  confirming  that  the  a-e- 
transition  does  not  occur.  The  velocity  and  pressure  profiles  are  quite  traditional  in  this 
case  and  agree  with  many  available  calculations  and  experiments  (e.g.,  with  relevant  ones 
published  in  [4]). 
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Next  calculations  shown  in  Figs.  2  and  3  correspond  to  Uo  =  0.5  km/s.  Results  shown  in 
Fig.  3  were  obtained  with  the  strength  data  for  e-phase  taken  from  Table  1.  The  shock 
profiles  in  these  calculations  evolve  till  T/  =  3.3  psec. 


0  x,  cm  2  0  x,  cm  2 

Figure  2:  Evolution  of  shock  profiles  at  Uo=0.5  km/s 

The  calculations  show  that  both  elasto-plastic  and  phase  transitions  take  place  in  this  case. 


0  X,  cm  2  0  X,  cm  2 


Figure  3:  Evolution  of  shock  profiles  at  Uo=0.5  km/s.  The  strength  of  e-phase  is  elevated. 

It  is  seen  that  the  strength  elevation  of  the  high-pressure  e-phase  does  influence  the 
profiles,  although,  not  very  significantly.  As  in  the  previous  calculation,  the  elasto-plastic 
transition  does  not  result  in  an  essential  temperature  change  and,  naturally,  gives  no  mass 
concentration  change.  The  mass  concentration  change  occurs  mainly  within  the  second 
transition  (a-e-transition).  Thus,  typical  three-wave  structure  is  observed  at  this  level  of 
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load.  It  is  interesting  to  note  that  the  profile  evolution,  when  developing  the  phase 
transition  within  the  shock  wave  with  time,  manifests  a  well-known  attenuation  of  the 
particle  velocity  and  pressure  behind  the  phase-transformation  step.  Discussion  on  this 
topic  can  be  found,  for  example,  in  [20]. 


Figure  4:  Evolution  of  shock  profiles  at  Uo=l  km/s 

Results  of  the  last  calculation  of  this  series  at  Lfo  =  1  km/ s  is  shown  in  Fig.  4  (T/  =  3.2  psec). 
In  this  case,  the  wave  velocity  of  the  phase-transforming  step,  while  increasing  with  Do, 
exceeds  the  wave  velocity  of  the  preceding  steps  and  the  three-wave  structure  degenerates 
into  a  single-wave  shock  wave  within  which  both  elasto-plastic  and  a-e-transitions  take 
place. 


Figure  5:  Free  surface  velocity  profiles  at  high  velocity  impact  by  Aluminium  plate 
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Last  calculations  aim  at  validation  of  the  model.  First  calculation  set-up  corresponds  to  an 
experiment  conducted  in  [21]  for  high  velocity  impact  of  2.46  mm-target  made  from  Iron 
by  2  mm  Aluminium  flyer  plate.  The  impact  velocity  is  1.9  km/ s.  Calculation  of  behaviour 
of  the  Aluminium  plate  is  conducted  with  the  model  and  data  taken  from  [4].  The  profiles 
drawn  in  Fig.  5  correspond  to  the  free  surface  velocity  of  the  target  taken  with  use  of  Laser 
Interferometry  (the  VISAR  technique).  Profile  E  in  Fig.  5  is  an  experimental  one  taken  from 
[21],  profile  Cl  is  calculated  with  the  EOS  and  constitutive  relation  data  used  in  the  test 
calculations  above,  and  profile  C2  is  calculated  with  the  elevated  strength  similarly  to  the 
calculation  in  Fig.  3.  An  extended  forerunner  in  front  of  the  elastic  precursor  of  the  profile 
E  is  explained  in  [21]  by  the  influence  of  an  attached  air  shock  being  developed  in  the  test. 
It  is  seen  from  Fig.  5  that  the  maximum  velocity  that  is  mainly  affected  by  properties  of  the 
e-phase  is  slightly  changed  when  comparing  Cl  and  C2,  but  the  change  is  not  dramatic 
(although  C2  agrees  better  with  the  experiment  E). 


Figure  6:  Free  surface  velocity  profiles  at  high  velocity  symmetric  impact  by  Iron  plate 

The  last  calculation  is  conducted  for  several  impact  velocities  in  the  symmetric  collision 
set-up  between  2  mm  iron  flyer  plate  and  4mm  iron  target.  The  target's  free  surface 
velocities  are  recorded  [22]  in  the  same  fashion  as  in  the  test  [21].  The  impact  velocities  are: 
1)  437  m/s;  2)  779  m/s;  3)  1079  m/s;  and  4)  1457  m/s.  The  corresponding  experimental 
profiles  (E1-E4)  taken  from  [22]  and  calculated  profiles  (C1-C4)  are  shown  in  Fig.  6.  The 
calculation  shows  a  good  agreement  and  a  pleasing  description  of  the  kinetic  features 
associated  with  the  wave  degeneration  when  impact  velocity  is  increasing. 
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9.  Conclusions 


The  present  work  has  generalised  the  decoupled  strength  model  [4]  for  the  case  of  two- 
phase  non-dispersive  materials.  In  doing  so,  the  decoupling  principle  was  preserved  that 
will  allow  implementation  of  the  model  in  commercial  hydrocodes.  Thermodynamic 
correctness  of  the  model  has  been  established  and  conditions  of  the  correctness 
formulated. 

Relations  for  the  constitutive  equations  of  the  model  have  been  built  up  for  Iron  and  a 
kinetic  of  the  bcc-to-hcp  phase  transition  of  Iron  described.  The  model  has  been  validated, 
using  experimental  data  on  the  a-c-phase  transition  available  in  literature.  Kinetic 
behaviour  of  the  elasto-plastic  and  phase  transition  resulting  in  a  three-wave  structure 
corresponding  to  both  transitions  has  been  observed  in  the  calculations  using  the  model. 
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